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Thermal dark matter that couples more strongly to electrons and photons than to neutrinos will 
heat the electron-photon plasma relative to the neutrino background if it becomes nonrelativistic 
after the neutrinos decouple from the thermal background. This results in a reduction in N e ff 
below the standard-model value, a result strongly disfavored by current CMB observations. Taking 
conservative lower bounds on N e ff and on the decoupling temperature of the neutrinos, we derive 
a bound on the dark matter particle mass of m x > 3 — 9 MeV, depending on the spin and statistics 
of the particle. For p-wave annihilation, our limit on the dark matter particle mass is stronger than 
the limit derived from distortions to the CMB fluctuation spectrum produced by annihilations near 
the epoch of recombination. 



Roughly 20—25% of the total energy content of the uni- 
verse is in the form of non-baryonic dark matter. While 
a dark matter particle mass in the GeV range is often as- 
sumed, there has also been interest in masses in the MeV 
range. Dark matter with a mass in this range was invoked 
to explain the 511 keV 7-rays observed by INTEGRAL 
[l| , and to explain the cosmic 7-ray background at 1 — 20 
MeV (2}- MeV dark matter can have interesting effects 
on large-scale structure [3| . Supersymmetric models with 
MeV dark matter have been proposed Q , and MeV dark 
matter can arise in the context of the WIMPless dark 
matter model [5J. 

We note here that a thermal MeV dark matter particle 
that couples more strongly to electrons and photons than 
to neutrinos will heat the electron-photon plasma when 
it becomes nonrelativistic before its abundance freezes 
out. If this occurs after the neutrinos decouple from 
the thermal background, then the ratio of the neutrino 
temperature to the photon temperature will be reduced, 
a process similar to the heating that occurs when the 
electron-positron pairs become nonrelativistic. The final 
result is a decrease in the effective number of neutrino de- 
grees of freedom. This effect was first explored by Kolb 
ct al. Q and more recently by Serpico and Raffelt Q in 
the context ofprimordial nucleosynthesis. Recent CMB 
observations |8l-flfj| place severe lower bounds on N e ff, 
allowing us to constrain this process. 

At recombination, the energy density in relativis- 
tic particles includes photons, whose temperature, T 7 , 
and therefore energy density is extremely well-measured, 
and a neutrino background with temperature T v — 
(4/ll) 1 / 3 T 7 . The theoretical prediction for the effective 
number of neutrinos (assuming slight reheating of the 
neutrinos from early e + e _ annihilation) is N e ff = 3.046 
[TTI . [l2| . The neutrino density cannot be measured di- 
rectly, but it can be inferred from measurements of the 
CMB. (For a discussion of the effect of N eff on the CMB 
fluctuations, see Refs. [H,[l3|)- The values of N e ff from 
recent CMB observations, in combination with other cos- 
mological data, are N ef} = 4.34±g"|| (68% CL) from 
WMAP @, N eff = 4.56 ± 0.75 (68% CL) from the At- 
acama Cosmology Telescope [§], and N e ff = 3.86 ± 0.42 
(68% CL) from the South Pole Telescope [Io|. Archidi- 
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ues of N e ff larger than the standard- model theoretical 
prediction, rather than smaller. 

The extent of the heating from dark matter annihila- 
tion in the early universe can be derived from entropy 
conservation (see Refs. [lfidi]]) from which our discus- 
sion is derived). Our paper assumes a dark matter par- 
ticle that couples much more strongly to electrons and 
photons than to neutrinos. The prototypical example of 
such a particle is one that interacts with ordinary mat- 
ter through an electromagnetic form factor, such as an 
electric or magnetic dipole (13 - 131) . 

Let XX denote the pair of dark matter particles. To 
make our study general, we will allow a range of possibili- 
ties for the dark matter, including a self-conjugate scalar, 
a non-self-conjugate scalar, a spin- 1/2 Majorana fermion 
or a spin-1/2 Dirac fermion. Thus, for the cases with 
self-conjugate and non-self-conjugate scalars, the nota- 
tion XX really means XX and XX* respectively. But for 
simplicity, we will keep the notation XX throughout the 
paper. 

Consider first the case where the dark matter annihi- 
lates entirely after the neutrinos decouple, which occurs 
at a temperature of T d « 2 - 3 MeV [ll|, |32|. The total 
entropy prior to XX annihilation is proportional to 
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For a relativistic particle, p = p/3, so following Ref. (l7j . 
we can write the total entropy density as 
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where <?*s is the total number of spin degrees of freedom 
for bosons, and 7/8 times the total number of spin de- 
grees of freedom for fermions. Then the total entropy 
is 

„ 2n 2 
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which is conserved through the process of any particle 
becoming nonrelativistic and annihilating. So the ratio 
of the final value of RT after annihilation to the initial 
value of RT prior to annihilation is 



(RT)f 
{RT)i 



g*si 
9*Sf 



1/3 
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where g*si and 3*s/ are the values of g*s for the rela- 
tivistic particles in thermal equilibrium before and after 
annihilation, respectively. When the XX pairs annihilate 
after neutrino decoupling, the neutrinos do not share in 
the heating, so that RT V is constant and T v oc while 
the photons and electron-positron pairs are heated as in 
Eq. ([5]). Therefore, for the XX pairs with g internal 
degrees of freedom, the ratio of T v to T 7 after XX anni- 
hilation is: 



if X is a fermion, and 
T„/T 7 = 
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if it is a boson. Taking, for example, the x particle to 
be a spin- 1/2 Major ana fermion gives g — 2, so that 
T„/T 7 = (22/29) 1/3 . Subsequent e+e" annihilation fur- 
ther heats the photon temperature relative to the neu- 
trino temperature by a factor of (11/4) 1 / 3 , so that the 



final ratio of the neutrino temperature to the photon tem- 
perature would be (88/319) 1 / 3 . 

In terms of N e ff, the energy density for neutrinos is 
given by 
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Since p v at fixed T 7 is the quantity that is inferred from 
CMB observations, a change in T u jT 1 will be interpreted 
as a change in N e ff, with N e ff oc (T u /T j ) 4l . In this case, 
XX annihilation reduces the value of T u /Ty relative to its 
value in the standard model by a factor of (22/29) 1 ^ 3 , 
which corresponds to N e ft = 3(22/29) 4 / 3 = 2.1, a value 
clearly excluded by the CMB observations. 

This value of N e f f corresponds to a dark matter parti- 
cle with a mass well below the neutrino decoupling tem- 
perature. However, to derive a useful limit, we must con- 
sider what happens when x annihilates during neutrino 
decoupling. Neutrino decoupling is not a sudden process, 
but for the purposes of our simplified calculation, we will 
take it to occur abruptly at a fixed temperature Td, and 
we will assume that dark matter annihilations before Td 
fully heat the neutrinos, while those after Td heat only 
the photons and e + e~ pairs. Let I(T 7 ) be given by (see, 
e.g., Weinberg 1972 for a similar calculation): 
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where the plus (minus) sign is for a fermionic (bosonic) 
dark matter particle, and the variable of integration is 
x = p x /T~. In the limit where all particles are fully 
rclativistic, I reduces to (27r 2 /45)3*s; the integral in Eq. 
([9]) just quantifies the contribution to / from XX as they 
become nonrelativistic. 

As mentioned above, the XX annihilation will heat up 
photons relative to neutrinos only after neutrino decou- 
pling. But this heating ends when the XX particles drop 
out of thermal equilibrium. Thus, the ratio of the neu- 
trino temperature to the photon temperature due to XX 
annihilation alone is 

1/3 
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where Tf is the temperature at which the XX particles 



freeze out. Since m x /Tf ~ 20 [i7j . it is obvious from Eq. 
© that we can simply set Tf = with negligible error: 



T,/T 7 
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The physical reason for this is that the XX abundance 
freezes out at a temperature of Tf ~ m x /20, while most 
of entropy from the XX annihilations is transferred to 
the thermal background when T ~ m x /3. Of course, the 
temperature ratio given by Eq. (|11[) must then be mul- 
tiplied by an additional factor of (4/11) 1 / 3 from e + e~ 
annihilations to obtain the final ratio of the neutrino tem- 
perature to the photon temperature. 

In this approximation, the effective number of neutri- 
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FIG. 1: The effective number of neutrino degrees of freedom, 
N e ff, that would be deduced from cosmic microwave back- 
ground observations for a thermal dark matter particle with 
mass m x , assuming sudden decoupling of the cosmic neutri- 
nos at a temperature Td- Curves correspond, top to bottom, 
to a g = 1 boson (short dash), g = 2 boson (dotted), g = 2 
fermion (solid) and g = 4 fermion (long dash). 



nos as measured by CMB experiments will be given by 



N eff = 3.046 



1(0) 
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I(T d ) 



(12) 



The value of N e ff as a function of m x /Td is shown in 
Fig. 1, for a self-conjugate scalar boson (g — 1), a non- 
self-conjugate scalar boson (g = 2), a spin-1/2 Majorana 
fermion (g = 2) and a spin-1/2 Dirac fermion (gf = 4). 

In fact, from Eqs. ©-(JT]), we can derive the m x -c Td 
limit for N e ff, namely 



N e ff = 3.046 
for fermionic \, and 

Neff = 3.046 
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for bosonic \- 

As noted earlier, neutrino decoupling is not a sudden 
process, so Td is not completely well-defined. Ref. [32[ 
gives a widely cited value of Td — 2.3 MeV for the elec- 
tron neutrinos, with the fi and r neutrinos decoupling at 



a higher temperature. However, neutrino oscillations will 
tend to equilibrate the decoupling of all three neutrinos, 
an effect discussed in Refs. [12|, |33[. Here we will simply 
take Td ^ 2 MeV as a conservative lower bound. 

Now we must determine a reasonable lower bound on 
N e ff. The combined results from Refs. [8l-[Io| are barely 
consistent with the standard model value of N e ff = 
3.046. However, we will err on the side of caution and 
choose a lower bound of N e ff > 2.6, which is excluded 
at 2cr by all three sets of CMB observations. 

These limits on N e f / and T d can be combined with the 
results displayed in Fig. 1 to derive a lower bound on m x . 
These bounds are m x ^ 3 MeV for the self-conjugate 
scalar boson, m x }t 6 MeV for a two-component boson 
or fermion, and m x ^ 9 MeV for a Dirac fermion. 

These limits are relevant for several models in the lit- 
erature. As noted by Beacom and Yuksel [34[ , the model 
proposed in Ref. |l[ actually requires positron injection 
at very low energies ( $, 3 MeV) toproduce the 511 keV 
7-rays observed by INTEGRAL [lj. But dark matter 
masses low enough to produce such particles from anni- 
hilations are ruled out by our limit. Thermal dark matter 
with the correct relic abundance interacting through an 
electric or magnetic dipole moment must have a mass less 
than 1 — 10 GeV to avoid conflict with direct detection 
experiments (28j ; our results shrink the allowed window 
from the other direction. 

Our limits are complementary to several others in the 
literature. As noted, dark matter particles with masses 
in this range also affect primordial nucleosynthesis, and 
bounds can be placed from the observed element abun- 
dances, particularly helium-4. However, the effect on 
N e f f as measured by the CMB appears to provide a bet- 
ter limit. For example, in the 1 — 10 MeV mass range, Ser- 
pico and Raffelt [3] found a maximum reduction of only 
0.002 in the primordial helium mass fraction. Using the 
results of Ref. [35[, this corresponds to AN e ff = —0.15, 
much smaller than the typical values in Fig. 1. How- 
ever, there is no contradiction between our results and 
those of Ref. 0. When T„/T 7 is reduced prior to pri- 
mordial nucleosynthesis, there are actually two effects on 
the helium-4 abundance. First, the reduction in the ex- 
pansion rate at fixed T 7 reduces the helium-4 abundance, 
and this is the dominant effect, as noted by Serpico and 
Raffelt. However, there is a second effect which partially 
cancels the first: the decrease in the electron neutrino 
temperature reduces the weak interaction rates, which 
tends to increase the hclium-4 abundance. Thus, the ef- 
fect on BBN is smaller than if one reduced the overall 
expansion rate alone. 

Another lower bound on m x comes from distortions to 
the CMB fluctuation spectrum due to annihilations near 
the epoch of recombination [36144 ll | . This effect excludes 
dark matter with masses ,$,1 — 10 GeV, a much tighter 
bound than ours (note that such annihilations also distort 
the spectrum of the CMB [42|, |43| , but these bounds are 
weaker given present observations). However, the CMB 
fluctuation bound only applies to s-wave annihilations, 
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for which (<jv) does not change between the dark matter 
particle freeze-out and the epoch of recombination. For 
p-wave annihilations, the annihilation rate at recombina- 
tion is generally negligible, and the CMB cannot be used 
to constrain such models. Our bounds, in contrast, do 
not depend on the velocity dependence of the annihila- 
tion cross section and therefore provide a good constraint 
in the case of p-wave annihilations. Indeed, the values of 
N eff derived in Refs. @-G3 assume a standard recombi- 
nation history, undistorted by dark matter annihilation, 
so it is unclear how s-wave annihilation at the epoch of 
recombination would affect the estimated values of N e f f . 
Of course, the reverse is also true; the bounds derived in 
Refs. [36 41] do not take into account the effect we have 
outlined in this paper. 

The bounds presented here can be evaded if the dark 
matter is asymmetric (see, e.g., Ref. [44| and references 
therein). Also, our bounds will be weakened to the extent 
that the dark matter couples to both the electron-photon 
plasma and to neutrinos. In fact, in the extreme oppo- 
site limit (coupling to neutrinos only), the \X annihila- 
tion heats the neutrinos instead of the photons, increas- 
ing N e ff and providing better agreement with current 
observations I4a|. 
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There is one obvious caveat to the bounds we have de- 
rived here. As noted earlier, the CMB limits on N e ft 
are only in marginal agreement even with the standard 
model value for N e ff. If future observations show con- 
clusive evidence that the observed N e ff disagrees with 
the standard model, some mechanism will be required to 
generate the additional rclativistic degrees of freedom, 
and this mechanism could also be invoked to erase the 
effects of the annihilating dark matter particle. Future 
PLANCK observations should help to resolve this issue. 
More precise observational bounds on N e ff would also 
justify a more exact treatment of the effect outlined here, 
going beyond our simplifying assumption of sudden neu- 
trino decoupling to a full numerical integration of the 
equations governing neutrino evolution in the early uni- 
verse. 
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